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Abstract 

We scan the parameter space of the NMSSM for the observability of at least one 
Higgs boson at the LHC with 300 fb~^ integrated luminosity, taking the present 
LEP2 constraints into account. We restrict the scan to those regions of parameter 
space for which Higgs boson decays to other Higgs bosons and/or supersymmetric 
particles are kinematically forbidden. We find that if lypl^-fusion detection modes 
for a light Higgs boson are not taken into account, then there are still significant 
regions in the scanned portion of the NMSSM parameter space where no Higgs bo- 
son can be observed at the 5ct level, despite the recent improvements in ATLAS 
and CMS procedures and techniques and even if we combine all non-fusion discov- 
ery channels. However, if the l^iy-fusion detection modes are included using the 
current theoretical study estimates, then we find that for all scanned points at least 
one of the NMSSM Higgs bosons will be detected. If the estimated 300 fb~^ sig- 
nificances for ATLAS and CMS are combined, one can also achieve 5a signals after 
combining just the non-VFVF-fusion channels signals. We present the parameters 
of several particularly difficult points, and discuss the complementary roles played 
by diff'erent modes. We conclude that the LHC will discover at least one NMSSM 
Higgs boson unless there are large branching ratios for decays to SUSY particles 
and/or to other Higgs bosons. 
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1 Introduction 



Supersymmetric extensions of the standard model generally predict relatively light Higgs 
bosons. One of the most important tasks of the LHC is the search for Higgs bosons 0. 
An important milestone in understanding the potential of the LHC was the demonstration 
that at least one Higgs boson of the minimal supersymmetric standard model (MSSM) 
would be detectable at the > 5a level throughout all of the MSSM parameter space so 
long as top squark masses do not exceed 1.5 to 2 TeV and so long as large branching 
fractions to decay channels containing supersymmetric particles are not substantial. 

In the present paper, we study, subject to these same and a few other simplifying 
restrictions, the detectability of Higgs bosons in the next-to-minimal supersymmetric 
standard model (NMSSM). In the NMSSM, one Higgs singlet superfield, S, is added to 
the MSSM in order to render unnecessary the bilinear superpotential term ^HiH2 by 
replacing it with \SH1H2, where the vacuum expectation value of the scalar component 
of S, (S), results in an effective bilinear Higgs mixing with fi = \{S). The detectability of 
the NMSSM Higgs bosons was first considered in a contribution to Snowmass 96 0]. The 
result, using the experimentally established modes and sensitivities available at the time, 
was that substantial regions in the parameter space of the NMSSM were found where 
none of the Higgs bosons would have been observable either at LEP2 or at the LHC even 
with an integrated luminosity of 600 fb^^ (two detectors with L = 300 fb~^ each). 

Since then, progress has been made both on the theoretical and the experimental sides. 
On the theoretical side, the dominant two-loop corrections to the effective potential of 
the model have been computed [|, |^ . These lead to a modest decrease in the mass of the 
lightest Higgs scalar, holding fixed the stop sector parameters. Inclusion of the two-loop 
corrections thus increases somewhat the part of the NMSSM parameter space excluded 
by LEP2 (and accessible at the Tevatron) but is of less relevance for the LHC. On 
the experimental side the expected statistical significances have been improved since 1996 
[|I], 3. Most notably, associated tih production with h ^ bb (originally discussed in 0), 
which in the SM context is particularly sensitive to ruh < 120 GeV, has been added by 
ATLAS and CMS to the list of Higgs boson detection modes [l|, ||, 0, . Analysis of this 
mode was recently extended [^] to ruh = 140 GeV, which, though not relevant in the SM 
case due to the decline in the bb branching ratio as the WW* mode increases, is highly 
relevant for points in our searches for which the WW* mode is suppressed in comparison 
to the SM prediction. In addition, techniques have been proposed for isolating signals 
for WW fusion to a light Higgs boson which decays to rr or WW^*\ 

It turns out that adding in just the tth process renders the no-Higgs-discovery param- 
eter choices described and plotted in including the "black point" described in detail 
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there, visible |jTT|]. In the present paper, we search for any remaining parameter choices 
for which no Higgs boson would produce a > So" signal. In this search, we perform a 
scan over nearly all of the parameter space of the model, the only parameter choices not 
included being those for which there is sensitivity to the highly model-dependent decays 
of Higgs bosons to other Higgs bosons and/or superparticles. The outcome is that, for 
an integrated luminosity of 300 fb~^ at the LHC, there are still regions in the parameter 
space with < 5a expected statistical significance (computed as Nsd = S/\fB for a given 
mode) for all Higgs detection modes so far studied in detail by ATLAS and CMS, i.e. 
including the tth tthb mode but not the lyW^-fusion modes. On the other hand, the 
expected statistical significance for at least one of these detection modes is always above 
3.6(T at 300 fb~^, and the statistical significance obtained by combining (using the naive 
Gaussian procedure) all the non-VTiy-fusion modes is at least 4.8cr. However, we find 
that all such cases are quite observable (at > 10. la) in one of the VrVT-fusion modes 
(using theoretically estimated statistical significances for these modes). For all points in 
the scan of parameter space, statistical significances obtained by combining all modes, 
including VTiy-fusion modes, are always >10.7cr. Thus, NMSSM Higgs discovery by just 
one detector with L = 300 fb~^ is essentially guaranteed for those portions of parame- 
ter space for which Higgs decays to other Higgs bosons or supersymmetric particles are 
kinematically forbidden. This represents substantial progress towards guaranteeing LHC 
discovery of at least one of the NMSSM Higgs bosons. 

In order to clarify the nature of the most difficult points in those portions of parameter 
space considered, we present, in sect. 4, examples of particularly difficult bench mark 
points for the Higgs sector of the NMSSM. Apart from the "bare" parameters of the 
model, we give the masses and couplings of all Higgs scalars, their production rates 
and branching ratios to various channels (relative to the SM Higgs) and details of the 
statistical significances predicted for each Higgs boson in each channel. The latter will 
allow an assessment of exactly what level of improvement in statistical significance will be 
required in the various different detection modes in order to render marginal modes visible. 
Of course, our estimates of the expected statistical significances are often somewhat crude 
(e.g. their dependence on the accumulated integrated luminosity). We beheve that our 
procedures always err in the conservative direction, leading to statistical significances that 
might be a bit small. Thus, the LHC procedures for isolating Higgs boson signals could 
provide even more robust signals for NMSSM Higgs boson detection than we estimate 
here. 

The detection modes, which serve for the searches for standard model or MSSM Higgs 
bosons, include (using the notation /i, a for CP-even, CP-odd Higgs bosons, respectively): 
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1) gg -fr, 

2) associated Wh or tih production with 77^^ in the final state; 

3) associated tth production with h — > bb; 

A) gg —>■ h/ a or associated bbh/a production with h/a ^ tt; 
5) gg ZZ'-*^ -> 4 leptons; 

Q) gg^h^ WW^*'^ l+l-yV- 

7) LEP2 e+e- Zh and e+e" ha; 

8) WW ^h^TT- 

9) WW WW^*\ 

where 8) and 9) are those analyzed at the theoretical level in and included in the 
NMSSM analysis for the first time in this paper. The above detection modes do not employ 
the possibly important decay channels i) h ^ hh, ii) h aa, iii) h h^h^ , iv) h — > aZ, 
v) a — > ha, vi) a —>■ hZ, vii) h,a h^W^, viii) h,a —>■ tt and ix) t —>■ h^b. The 
decay modes i)-vii) give high multiplicity final states and deserve a dedicated study [|12 



while the existing analyses of the tt final state signatures are not very detailed. Further, 
when kinematically allowed, the t h^b signal would be easily observed according to 
existing analyzes. Thus, in this paper we restrict our scan over NMSSM parameter space 
to those parameters for which none of these decays are present. In addition, we take the 
constraints of LEP2 [via the mode 7)] into account, and only accept points for which 5a 
discovery at LEP2 would not have been possible |1^ . 

The Higgs sector of the NMSSM consists of 3 scalars, denoted /ii, /i2, ^3 with < 
rrih^ < rrih^, 2 pseudo-scalars, denoted 01,02 with m^^ < m^j, and a charged Higgs pair, 
denoted h^. Mixing of the neutral doublet fields with the gauge singlet fields in the 
scalar and in the pseudo-scalar sector can be strong. The scalar mixing can lead to 
a simultaneous suppression of the couplings of all the hi to gauge bosons, and hence 
to a suppression of many of the detection modes above. (Of course, the have no 
tree-level couplings to gauge boson pairs and the one-loop couplings are too small to 
yield useful event rates.) The couplings of the Higgs bosons to t- or b-quarks can be 
amplified, reduced or even change sign with respect to the standard model couplings. 
Hence negative interferences can occur among the (loop-) diagrams contributing to gg 
hi and hi — > 77, leading again to suppressions of the above detection modes. A complete 
simultaneous annulation of all detection modes is not possible, but simultaneous reduction 
of all detection modes is possible and it is for such parameter choices that NMSSM Higgs 
boson discovery is most difficult. 

In the next section, we define the class of models we are going to consider, and the way 
we perform the scan over the corresponding parameter space. In section 3 we describe 
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our computations of the expected statistical significances of the detection modes 1) - 
9) above. In section 4, we present six particularly difficult bench mark points (in table 
1) and details regarding their statistical significances in channels l)-9) in table 2, with 
a summary of overall statistical significances in table 3. Using these tables, we give a 
discussion of the properties of these points. 

2 NMSSM Parameters and Scanning Procedure 



In this paper, we consider the simplest version of the NMSSM |T^, 0, ^7j, where the term 
HH1H2 in the superpotential of the MSSM is replaced by (we use the notation A for the 
superfield and A for its scalar component field) 

XH1H2S + , (2.1) 

so that the superpotential is scale invariant. We make no assumption on "universal" soft 
terms. Hence, the five soft supersymmetry breaking terms 

ml^Hf + ml^Hl + m^^^ + XAxH^H2S + ^A^S^ (2.2) 

are considered as independent. The masses and/or couplings of sparticles are assumed 
to be such that their contributions to the loop diagrams inducing Higgs production by 
gluon fusion and Higgs decay into 77 are negligible. In the stop sector, which appears 
in the radiative corrections to the Higgs potential, we chose the soft masses mq = ttit = 
Msusy = 1 TeV, and varied the stop mixing parameter 

X = 2 II — I (2 3) 

Mi,, + ml V 12(Mi,, + m?) y • ^ ' ^ 

As in the MSSM, the value Xt = - so called maximal mixing - maximizes the radiative 
corrections to the Higgs masses, and we found that it leads to the most challenging points 
in the parameter space of the NMSSM. 

Assuming that the Higgs sector is CP conserving, the independent parameters of the 
model are thus: A, k, m'jj_^, mjj^, m|, A\ and A^. For purposes of scanning and analysis, it 
is more convenient to eliminate mj^^, mj^^ and m| in favor of Mz, tan/5 and fics = X{S) 
through the three minimization equations of the Higgs potential (including the dominant 
1- and 2- loop corrections [§) and to scan over the six independent parameters 

A, K, tan P, /ieflf, ^A, • (2.4) 

We adopt the convention A, k > 0, in which tan/5 can have either sign. The absence of 
Landau singularities for A and k, below the GUT scale (~ 2 x 10^^ GeV) imposes upper 
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bounds on these couplings at the weak scale, which depend on the value of ht and hence 
of tan/? |jT5|, |16|. Using m^^'^ = 175 GeV, one finds A 

max ~ 0.69 and /tmax ~ 0.62 for 



intermediate values of tan (3. 

For each point in the parameter space, we diagonalize the scalar and pseudo-scalar 
mass matrices and compute the scalar, pseudo-scalar and charged Higgs masses and cou- 
plings taking into account the dominant 1- and 2- loop radiative corrections ||^. We then 
demand that the Higgs scalars satisfy the LEP2 constraints on the e+e~ —>■ Zhi produc- 
tion mode (taken from |]13|, fig. 10), which gives a lower bound on m/j^ as a function of 



the ZZhi reduced coupling. We also impose LEP2 constraints on e~^e hiaj associated 
production (from |T^, fig. 6), yielding a lower bound on m/j. + as a function of the 



Zhittj reduced coupling. 

In order to render the above-mentioned processes i) - ix) kinematically impossible, we 
require the following inequalities among the masses: 

m,j3 < 2m h^, 2ma^, 2m h±, ma^ + Mz, mh± + Mw] 
ma^ < mh^ + ma^, m^^ + M^, m^i + Mw] > 155GeV. 

In addition we require |/ieff| > 100 GeV; otherwise a light chargino would have been 
detected at LEP2. (The precise lower bound on \fies\ depends somewhat on tan/5 and 
the precise experimental lower bound on the chargino mass; however, our subsequent 
results do not depend on the precise choice of the lower bound on |/iefr|-) We further 
note that for the most challenging parameter space points that we shall shortly discuss, 
\Hes\ > 100 GeV is already sufficient to guarantee that the NMSSM Higgs bosons cannot 
decay to chargino pairs so long as the SU(2) soft-SUSY-breaking parameter M2 is also 
large. In fact, in order to avoid significant corrections to 77/14 and 7704 couplings coming 
from chargino loops it is easiest to take M2 ^ fies (or vice versa). This is because the 
hixtXi coupling is suppressed if the xt is either pure higgsino or pure gaugino. Since 
the parts of parameter space that are challenging with regard to Higgs detection typically 
have ~ 100 — 200 GeV, the validity of our assumptions requires that M2 be large and 
that the chargino be essentially pure higgsino. 

Using a very rough sampling, we determined, as expected from previous work 0, that 
it is only for moderate values of tan /3 that < So" signals might possibly occur. From this 
sampling, we determined the most difficult parameter space regions and further refined 
our scan to the following: 

• 4.5 < I tan/5| < 8 (both signs) in steps of 0.25; 

• 0.001 < A < min[0.21, Amax], using 20 points; 
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• 0.001 < K < min[0.24, Kmax], using 20 points; 

• 100 GeV < l/ioffl < 300 GeV (both signs), in steps of 5 GeV; 

• < < 160 GeV, with Ax opposite in sign to /igfr, using steps of 5 GeV; 

• 25 GeV < \Af^\ < 170 GeV, with opposite in sign to /lefr, using steps of 5 GeV. 

For those points sampled in this final scan which satisfy all the constraints detailed earlier, 
we compute the expected statistical significances for the processes 1) to 9) listed in section 
1, as described in the next section. As a rough guide, from the ~ 10^ points detailed in 
the above list, we find about 250,000 that pass all constraints and have Nsd < 5 (for 
L = 300 fb~^) in each of the individual discovery modes 1) - 7). We shall tabulate a 
number of representative points taken from this final set in section 4. 



3 Expected Statistical Significances 

From the known couplings of the NMSSM Higgs scalars to gauge bosons and fermions 
it is straightforward to compute their production rates in gluon-gluon fusion and various 
associated production processes, as well as their partial widths into 77, gauge bosons and 
fermions, either relative to a standard model Higgs scalar or relative to the MSSM H 
and/or A. This allows us to apply "NMSSM corrections" to the processes 1) - 9) above. 

These NMSSM corrections are computed in terms of the following ratios. For the 
scalar Higgs bosons, cy is the ratio of the coupling of the hi to vector bosons as compared 
to that of a SM Higgs boson (the coupling ratios for hiZZ and hiWW are the same), 
and Q, Cfc are the corresponding ratios of the couplings to top and bottom quarks (one 
has Cr = Cb). Note that we always have |cy| < 1, but q and Cf, can be larger, smaller or 
even differ in sign with respect to the standard model. For the CP-odd Higgs bosons, cy 
is not relevant since there is no tree-level coupling of the to the VV states; q and Cb 
are defined as the ratio of the 275 couplings for tt and bb, respectively, relative to SM-like 
strength. 

We emphasize that our procedure implicitely includes QCD corrections to the Higgs 
production processes at precisely the same level as the experimental collaborations. First, 
the ATLAS and CMS collaborations employed Monte Carlo programs such as ISAJET 



18| and PYTHIA |19| in obtaining results for the (MS)SM. These programs include 
many QCD corrections to Higgs production in a leading-log sense. This is the best that 
can currently be done to implement QCD corrections in the context of experimental 
cuts and neural-net analyses. Clearly the more exact NNLO results for many of the 
relevant processes will slowly be implemented in the Monte Carlo programs and increased 
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precision for Higgs discovery expectations will result. Since our goal is to obtain NMSSM 
results that are completely analogous to the currently available (MS)SM results, we have 
proceeded by simply rescaling the available (MS)SM experimental analyses. In doing the 
rescaling of the Higgs branching ratios we have included all relevant higher-order QCD 
corrections ^ using an adapted version of the FORTRAN code HDECAY pl|. f\ We 
now give additional details on our rescaling procedures. 

The expected statistical significances for the processes 1) and 6) are computed be- 
ginning with results for the SM Higgs boson taken from ref fig.l ("Expected Ob- 
servability of Standard Model Higgs in CMS with 100 fb~^"). The application of the 
NMSSM corrections using cy, Ct and Cb [which determine T{gg — > hi), BR{hi — > 77) and 
BR{hi WW*)] is straightforward in these two cases. 

The expected statistical significances for process 2) are taken from the same figure. 
In ref. |^ one finds that Whi and tihi production contribute with roughly equal weight 
to the SM signal. This allows us to decompose the expected significance into the corre- 
sponding production processes, apply the NMSSM corrections, and then recombine the 
production processes. 

The expected Standard Model Higgs statistical significances for process 3) are taken 
from table 19-8 in ref. [0], with the extension to Higgs masses above 120 GeV as pro- 
vided in 0, using a numerical interpolation for Higgs masses below 140 GeV. For the 
standard model process 5) we again use ref. |jl|, tables 19-18 and 19-21. In both cases, 
the application of the NMSSM corrections is straightforward. 

The estimation of the statistical significances for the process 4) in the NMSSM requires 

the most discussion. Figure 19-62 of ref. [jl[ and fig. 8 of ref. give the 5cr contours 

in the tan (3 - tua plane of the MSSM. The critical issue is how much of these So" signals 

derive from gg ^ H + gg ^ A production and how much from associated bbH + bbA 

production, and how each of the gg fusion and bb associated production processes are 

divided up between H and A. For the former, we turn to table 19-35 of ref. [l|]. There, 

we see that it is for cuts designed to single out the associated production processes that 

large statistical significance can be achieved and that such cuts provide 90% of the net 

statistical significance of Nsd = 8.9 (3.9 for gg fusion cuts combined in quadrature with 

8.0 for bbH+bbA associated production cuts) for = 150 GeV and L = 30 fb^^. (For the 

tin our computations however, we neglect the contribution to the hij'j coupHng coming from the 
charged Higgs loop. Despite the relatively small masses of the for our most problematical points, 
the charged Higgs loop decouples especially for the small values of A and XA\ characteristic of 
difficult points for which the actual h^h~hi coupling is only of order a few times gmvy 7(4^2) Its 
contribution would typically only be of order a few percent even though our difficult points have smaller 
77 /li coupling than a SM-like Higgs boson by virtue of suppressed hiWW coupling and/or cancellation 
between the top and W loop contributions. 
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associated production cuts, the table shows that the contribution of the gg fusion processes 
to the signal is very small.) The percentage of Nsd deriving from (yf^f-fusion cuts is even 
smaller at high mA- Since we are mainly interested in m/f, tua G [100 GeV, 200 GeV], 
we will assume that 90% of the statistical significance along the contours of fig. 19-62 
comes from the associated production cut analysis; this will give us a slightly conservative 
estimate of the associated production Nsd values at still higher uia- With this choice, 
the 5cr contour at L = 100 fb^^ from fig. 19-62 of ref. [jl| corresponds to a 4.5o" contour 
for associated bbH + bbA production alone. Since the values of tan P along this contour 
are large, we can separate the H and A signals from one another by using the following 
properties of the MSSM within which fig. 19-62 of ref. |I| was generated: (a) BR{H 
tt) ~ BR{A — s> rr) ~ 0.09; (b) the hbA and hbH couplings are very nearly equal and 
scale as tan/5; and (c) ~ mn within the rr mass resolution. As a result, the net 
signal rate along this contour is approximately twice that for hbA or hbH alone. Thus, 
Nsd = 2.25 would be achieved for hbA or hbH along this contour were uia and mu widely 
separated. 

We can then compute the statistical significance for the hbhi and hbai signals with 
hi, a-i TT decay using the following procedure. First, the NMSSM hbhi and hbai pro- 
duction rates are related to the MSSM hbH and hbA rates by the factors [cf,(/;.j)]^/ tan^ (3 
and [cb(aj)]^/ tan^/3, respectively. Next, we account for the fact that the rr branching 
ratios of the NMSSM scalars and pseudo-scalars differ somewhat from the value of 0.09 
appropriate for the MSSM H and A. In particular, BR{h^ tt) is significantly reduced 
when the h^ has large enough mass and large enough cy that it acquires a modest WW* 
branching ratio. Typical reductions will be tabulated in table 2. Thus, defining the value 
of tan/? as a function of shown by the 100 fb~^ curve of fig. 19-62 in ref. |1| as 
tan/52.25('TiA), we compute NsD{hi) for L = 100 fb^"*^ as 

1 2 



NsDihi) = 2.25 



Cb{hi 



BRMhi), (3.1) 



_tan/52.25("^h, 

where BR^t = BR{hi tt)/BR{H tt) = BR{hi -> rr)/0.09. An exactly parallel 
procedure is employed for NsD{(ii)- 

The above procedure is conservative in that it assumes no contribution to the rr 
channel Nsd from the gg fusion processes. However, as we shall describe in the next 
section, for the most difficult points in parameter space the gg-iusion rates are very 
substantially suppressed relative to MSSM values. For these points, essentially 99% of 
the rr channel Nsd would derive from 66-|-Higgs associated production. 

In recombining the scalar and pseudo-scalar signals, we must account for the fact that 
they can have fairly different masses in the NMSSM. In this paper, we have chosen to 
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recombine the scalar and pseudo-scalar signals at different masses following the procedure 
of ref. section 5.4, with am ~ 30 GeV as estimated from fig. 19-61 in [|l| at 



high luminosity and extrapolated to rriA ^ 150 GeV. This procedure leads to somewhat 
approximate estimates of the NMSSM statistical significances for this detection mode. 

Results for the statistical significances of the hi signals in modes 8) and 9) were simi- 
larly obtained by rescaling the theoretical results of [1^ ( summarized most conveniently 
in the last of the listed papers) using the values of [cv{hi)]'^, BRtt and BRWW*. 

Using the above procedures, for each point in the parameter space of the NMSSM we 
obtain the statistical significances predicted for an integrated luminosity of 100 fb~^ for 
each of the detection modes 1) - 9). In order to obtain the statistical significances for the 
various detection modes at 300 fb^^, we multiply the 100 fb~^ statistical significances by 
a/3 in the cases 1), 2), 3), 5) and 6), but only by a factor of 1.3 in the cases 4), 8) and 9). 
That such a factor is appropriate for mode 4), see, for example, fig. 19-62 in [|l|. Use of 
this same factor for modes 8) and 9) is simply a conservative guess. 



4 Difficult Points 

As stated in the introduction we still find "black spots" in the parameter space of the 
NMSSM, where the expected statistical significances for all Higgs detection modes 1) - 7) 
are below 5a at 300 fb~\ The reasons for this phenomenon have been described above; 
see also the corresponding discussion in 0. However, after including the modes 8) and 
9), the points that provide the worst 1) - 6) statistical significances typically yield robust 
signals in one or the other of the lyiy-fusion modes 8) and 9). 

In order to render the corresponding suppression mechanisms of the detection modes 
reproducible, we present the detailed properties of several difficult points in the parameter 
space in table 1. The notation is as follows: The bare parameters are as in eq. (2.5), 
with m^^, and m| fixed implicitly by the minimization conditions. (As noted earlier, 
with the convention A, k > in the NMSSM, the sign of tan (3 can no longer be defined 
to be positive.) For the reasons discussed below eq. (|2.3|) we chose in the stop sector 
niQ = rriT = Msusy = 1 TeV and Xf = VG for all of the points (1 - 6). We have also 
fixed mfop"^ = 175 GeV. For both scalar and pseudoscalar Higgs bosons, "gg Production 
Rate" denotes the ratio of the gluon-gluon production rate with respect to that obtained 
if Ct = Cfe = 1, keeping the Higgs mass fixed. For scalar hi, this is the same as the ratio 
of the gg production rate relative to that predicted for a SM Higgs boson of the same 
mass. For the scalar hi, BR'j'j denotes the ratio of the 77 branching ratio with respect to 
that of a SM Higgs boson with the same mass. (A verification of the reduced gluon-gluon 
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production rates or 77 branching ratios would sometimes require the knowledge of the 
couplings to higher precision than given, for convenience, in table 1.) Also given for the 
scalar hi are the ratios BRbb and BRWW* of the bb and WW* branching ratios relative 
to the SM prediction (as noted above, one has BRtt — BRbb). 

In table 2, we tabulate the statistical significances for the hi in all the channels 1) 
- 9); production of the CP-odd turns out to be relevant only when they add to the 
hi signals in process 4). Also note that, all these problematical points are such that 
m/jj + > 206 GeV, so that e'^e~ — > /ii + Oi followed by /ii, Oi — > bb would have 
been kinematically forbidden at the highest LEP2 energy. Hence, for LEP2 mode 7) we 
only give the statistical significance for e+e~ — > Zhi. Also tabulated in table 2 are four 
statistical significances obtained by combining various channels. This combination is done 
in the Gaussian approximation: 

1/2 

(4.1) 

where runs over the channels i being combined. We give results for the following 
combinations: 

a) NsD obtained by combining LHC channels 1) - 6); 

b) NsD obtained by combining LHC channels 1) - 6) and LEP2; 

c) ^SD obtained by combining LHC channels 1) - 6) with the H^lV-fusion channels 8) 

and 9); 

d) NsD obtained by combining all LHC channels and LEP2, i. e. by combining all channels 

l)-9). 

In those cases where there is no LEP2 signal, a)=b) and c)=d). In addition, in our point 
selection we have required a mass difference of at least 10 GeV between scalar Higgses, 
so that they yield well separated signals and no statistical significance combination of 
two different scalar Higgses is needed. All parameter choices for which Higgs boson 
masses differ by less than 10 GeV yield stronger signals than the cases retained. (The 
increased net signal strength of overlapping Higgs signals in those channels with limited 
mass resolution arises as a result of Nf^{l + 2) - {Si + S2)/^/B > yjsf + Si/y/B.) 

As summarized in table 3, all of the tabulated "bench mark points" have statistical 
significances below 5a for all of the detection modes 1) - 6) at 300 fb~^ and 7) at LEP2. 
In more detail, as tabulated in table 2 and summarized in table 3, the best signals in the 
modes 1) - 6) for the points #1 - #6 at the LHC are: 

• point ^1, NsD =4.37 for mode 2) and hi, 



Arcombined 
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• point #2, NsD =3.95 for mode 3) and h2] 

• point #3, NsD =3.62 for mode 4) and h^, 

• point #4, NsD =4.46 for mode 5) and h^, 

• point #5, NsD =4.83 for mode 3) and hi, 

• point #6, NsD =4.86 for mode 4) and h^; 

Further, for point #3, the combined statistical significance of modes 1) - 6) (also tabulated 
in table 3) would still be below 5 for any one hi, although v^N^-^^ > 5 (as is likely to 
be relevant by combining ATLAS and CMS data once each detector has accumulated 
L = 300 fb"^) for at least one of the h^. However, for all these "difficult" points the 
VrVT-fusion modes 8) and/or 9) provide (according to theoretical estimates) a decent 
(sometimes very strong) signal. 

The points #1 - ^^4 differ as to which of the modes 1) - 6) and which hi yields 
the largest statistical significance should the lyi^-fusion mode 8) not provide as strong 
a signal as suggested by the theoretical estimates. To render these points observable 
without the VTl^-fusion mode 8) would require improvements of all detection modes 2) 
-5). 

As in 10, we find that difficult points in the parameter space generally have | tan j3\ ~ 5. 
This is the region of tan (3 for which the bbh, bba signals are still not very much enhanced 
but yet the gg —>■ h,a and tth,tta signals have been suppressed somewhat. In a few 
cases, however, difficulties also arise for |tan/5| as large as 8, as shown in the case of 
point #5. Also as in ||^, the most difficult points are those in which the masses of the hi 
and are relatively close in magnitude, typically clustered in a ~ 60 GeV interval above 
~ 105 GeV. Such clustering maximizes the mixing among the different Higgs bosons 
and thereby minimizes the significance of the discovery channels for any one Higgs boson. 
In particular, it is for strong mixing among the hi that the statistical significance for 
discovery modes based on a large VV coupling for any one hi are most easily suppressed. 

Finally, for point #6, we have minimized the statistical significances for the WW- 
fusion modes over the parameter space, while keeping the statistical significances of modes 
1) - 6) below 5. One can see that it still gives a strong lO.lcr signal in mode 8). [Smaller 
NsD for mode 8) would have been possible if we had allowed stronger signals in modes 1) 
- 6), in particular had we allowed smaller mass separation, < 10 GeV, between the two 
lightest Higgs bosons.] In addition, for point #6 = 112 GeV and the ZZ coupling of 
hi is sufficiently large that it would have yielded a 4.8cr signal at LEP2. Had we taken 
a top quark mass slightly larger, m^^'^ = 178 GeV, we would have found a very similar 
point with a hi mass of ~ 115 GeV, which could have been responsible for the excess 



observed at LEP2 26 



12 



5 Discussion and Conclusions 



In this paper, we have addressed the question of whether or not it would be possible to fail 
to discover any of the Higgs bosons of the NMSSM using combined LEP2 and LHC data, 
possibly resulting in the erroneous conclusion that Higgs bosons with masses below 200 
GeV have been excluded. We have demonstrated that, assuming that the decay channels 
i) - ix) are either kinematically disallowed or render a Higgs boson observable, this is 
unlikely (at the > 5a level) to happen. Certainly, there are points in NMSSM parameter 
space for which the statistical significances for the individual detection modes 1) - 6) 
{i.e. those analyzed in detail by ATLAS and CMS) are all well below 5cr for integrated 
luminosity of 300 fb~^. However, by combining several of the modes 1) - 6) and 300 fb^^ 
data from both ATLAS and CMS, a > So" signal can be achieved based just on modes 1) - 
6). Further, we have found that throughout all of the NMSSM parameter space (scanned 
subject to the earlier listed restrictions) for which such weak signals in modes 1) - 6) 
are predicted, the theoretical estimates for the VTVT-fusion modes indicate that an easily 
detected WW —^h^rr signal should be present. Thus, our conclusion is that for all 
of the parameter space of the NMSSM compatible with reasonable boundary conditions 
for the parameters at the GUT scale (with, of course, non-universal soft terms in general) 
and such that Higgs pair and SUSY pair decays of the Higgs bosons are kinematically 
forbidden, at least one of the NMSSM Higgs bosons will be detected at the LHC. This 
is a big improvement over the results from the earlier Snowmass 1996 study which was 
somewhat negative without the inclusion of the tih —>■ tibb mode 3), and the W^VT-fusion 
modes 8) and 9). 

It is amusing to note that all of our bench mark points for which Higgs discovery 
is most difficult at the LHC include a Higgs scalar with mass close to 115 GeV (with, 
however, reduced couplings to the Z boson), which could be responsible for the excess 
observed at LEP2 pef. 

Another important point that appears from our analysis is the fact that the full L = 
300 fb~^ of integrated luminosity (per detector) is needed in order to have robust NMSSM 
Higgs discovery in the portion of parameter space considered here. Of course, as in the 
MSSM, it is very possible that only one of the CP-even NMSSM Higgs bosons might be 



detected at the LHC but that, as studied by Kamoshita et al. in |]T^, the observation of 
all the CP-even Higgs bosons of the NMSSM would be possible at the LC by virtue of all 
having some non-negligible level of ZZ coupling and not having very high masses. Even 
at the LC, the CP-odd Higgs bosons might escape discovery, although this would not be 
the case for the parameter choices that we have found which make LHC discovery of even 
one NMSSM Higgs bosons most challenging. This is because, for such parameters, the 
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ai are relatively light and could be readily seen at the LC in the processes e^e hiaj, 
e^e^ uVaitti and e+e" ^ Z* ^ Za^ai, assuming an integrated LC luminosity of 
1000 fb~^ and energy ^ > 500 GeV ^ . 



This study makes clear the importance of continuing to expand the sensitivity of 
existing modes and continuing to develop new modes for Higgs detection at the LHC in 
order not to have to wait for construction of a linear e+e~ collider for detection of at 
least one of the SUSY Higgs bosons. In particular, study of modes i) - ix) and SUSY 
pair channels should all be pushed. The problematical points that we have emphasized 
here are unlikely to be substantially influenced by ti or SUSY decays since all the Higgs 
masses are below ~ 200 GeV so that ti decays will be kinematically highly suppressed 
(one of the top quarks would have to be virtual) and SUSY pair decays are quite unlikely 
to be significant given LEP2 limits on the masses of SUSY particles. However, by allowing 
Higgs (in particular, pseudoscalar) masses such that one or more of the channels i)-vii) 
are kinematically allowed we have found points for which discovery in modes l)-9) will 
not be possible |l^. Thus, a full "no-lose" theorem for NMSSM Higgs boson discovery 
at the LHC will require exploring additional discovery modes sensitive to those portions 
of parameter space for which Higgs decays to other Higgs bosons are important, and 
might necessitate combining results from both the ATLAS and CMS detectors and/or 
accumulating more integrated luminosity. 
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Table Captions 

Table 1: We tabulate the input bare model parameters, the corresponding Higgs masses, 
and the corresponding Higgs couplings, relative to SM Higgs boson coupling strength, for 
6 bench mark points. Also given for the CP-even hi are ratios of the gg production rate 
and various branching fractions relative to the values found for a SM Higgs of the same 
mass. For the CP-odd Oj, "gg Production Rate" refers to the value relative to what would 
be found if both the 66 and the tt 75 couplings had SM-like strength. 

Table 2: Scalar Higgs statistical significances, Nsd — S/\/B, in various channels for the 
6 bench mark points. For each individual Higgs, we give (in order): Nsd for the channels 
1) - 9) described in the text; Gaussian combined Nsd for non-Vl^l^-fusion LHC channels; 
combined Nsd for non-VFVF-fusion LHC channels plus LEP2; combined Nsd for all LHC 
channels, including the fusion channels WW h ^ tt and WW ^ /i — > VFVF^*) 
channels; and combined Nsd for all LHC channels plus LEP2. 

Table 3: Summary for all Higgs bosons. The entries are: maximum non-l^VF fusion 
LHC Nsd; maximum LHC WW fusion Nsd', best combined Nsd after summing over 
all non-VFl^-fusion LHC channels; and best combined Nsd after summing over all LHC 
channels. The Higgs boson for which these best values are achieved is indicated in the 
parenthesis. One should refer to the preceding table in order to find which channel(s) 
give the best values. 
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Table 1 



Point Number 


1 


2 


3 


4 


5 


6 


Bare Parameters 




A 


0.0340 


0.0450 


0.0230 


0.0230 


0.1330 


0.0230 




0.0198 


0.0248 


0.0129 


0.0069 


0.1459 


0.0114 


tan/3 


6.00 


5.25 


-5.5 


5.75 


-8 


-6 




140 


-110 


115 


-235 


100 


150 


Ax{GeV) 


-35 


25 


-95 


40 


-135 


-100 


A,-[G( \ ) 


-130 


70 


-90 


80 


-75 


-110 


Scalar Masses and Couplings 




ruH, (GeV) 


115 


100 


103 


113 


114 


112 


cv 


-0.66 


0.32 


-0.34 


0.67 


-0.87 


-0.71 


Ct 


-0.65 


0.30 


-0.31 


0.65 


-0.81 


-0.66 


Cb 


-1.07 


0.66 


-1.27 


1.16 


-4.50 


-2.40 


gg Production Rate 


0.39 


0.08 


0.08 


0.39 


0.56 


0.36 




0.43 


0.26 


0.09 


0.38 


0.05 


0.11 


BRbb = BRtt 


1.12 


1.08 


1.10 


1.12 


1.18 


1.15 


BRWW^*'> 


0.42 


0.25 


0.08 


0.37 


0.04 


0.10 


TO/12 (GeV) 


125 


114 


114 


126 


144 


122 


Cv 


-0.74 


-0.83 


0.79 


-0.73 


0.46 


0.59 


Ct 


-0.72 


-0.74 


0.70 


-0.71 


0.57 


0.54 


Ch 


-1.49 


-3.28 


3.46 


-1.47 


-6.66 


2.24 


gg Pr(xlucti(jii Rate 


0.4(i 


0.44 


0.40 


0.4.-. 


1.18 


0.23 




0.33 


0.08 


0.07 


0.34 


0.01 


0.10 


rSiioO = rfKTT 


1.30 


1.18 


1.18 


1.32 


3.06 


1.31 


rSHW W ^ ' 


0.32 


0.08 


0.06 


0.33 


0.01 


0.09 


(GeV) 


205 


153 


148 


201 


202 


155 


cy 


-0.14 


-0.46 


-0.51 


-0.15 


0.18 


-0.39 


Ct 


-0.30 


-0.63 


-0.67 


-0.32 


0.17 


-0.55 


Cb 


5.80 


4.17 


4.20 


5.53 


0.68 


5.12 


gg Production Rate 


0.31 


0.84 


0.95 


0.33 


0.02 


0.80 


ID JX^^ 


0.13 


0.05 


0.05 


0.15 


0.98 


0.03 


JDnUU — JDJXTT 


308.66 


5.83 


3.92 


274.41 


13.97 


8.12 




0.18 


0.07 


0.06 


0.21 


0.96 


0.05 


Pseudo-Scalar Masses and Couplings 




ma, (GeV) 


191 


112 


130 


130 


113 


145 


Ct 


0.03 


-0.03 


-0.10 


-0.01 


-0.10 


-0.16 


Cb 


1.16 


-0.83 


-2.95 


-0.19 


-6.55 


-5.77 


gg Production Rate 


0.00 


0.00 


0.03 


0.00 


0.31 


0.08 


TOas (GeV) 


206 


141 


137 


198 


174 


158 




0.16 


0.19 


-0.15 


0.17 


-0.07 


-0.05 


Cb 


5.89 


5.18 


-4.64 


5.75 


-4.59 


-1.65 


gg Production Rate 


0.02 


0.07 


0.06 


0.02 


0.03 


0.00 


Charged Higgs Mass 




TOc(GeV) 


221 


162 


157 


213 


157 


167 
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Table 2 



Point 


1 


2 


3 


4 


5 


6 


(! n H "nnol 






"1 J^^ag 








Nsd{1) 


3.74 


0.35 


0.13 


3.18 


0.62 


0.83 




4.37 


0.59 


0.22 


3.92 


0.85 


1.22 


Nr. r^('i) 


9 7Q 


u.oo 


u.oo 


o.uo 




o.ou 


Nsd(^) 


0.08 


0.07 


0.76 


0.09 


4.52 


0.40 


Nsd{5) 


0.83 


0.00 


0.00 


0.64 


0.12 


0.16 


NsoiQ) 


1.10 


0.09 


0.03 


0.90 


0.16 


0.22 


Nsd{7) 


0.00 


3.37 


3.40 


3.29 


0.00 


4.79 


Nsd{8) 


9.29 


1.22 


1.59 


8.93 


16.78 


10.08 




2.39 


0.00 


0.00 


1.74 


0.41 


0.49 


\/EtiWsD{iW 

\l 1 J I — i L \ /J 


6.54 


1.09 


1.17 


5.99 


6.69 


3.65 


— ^7 Tat /'MO 

VE^=l[^SD(^)]2 


6.54 


3.55 


3.59 


6.84 


6.69 


6.02 




11.61 


1.64 


1.97 


10.89 


18.07 


10.73 




11. ()1 


3. 75 


3.93 


11.38 


18.07 


11.75 


I , n TiTi P 






"-2 nigg 








iVsi3(l) 


3.69 


0.83 


0.61 


3.62 


0.22 


0.55 




4.01 


1.25 


0.92 


3.93 


0.05 


0.74 






o.yo 


o.Oo 


Z.oU 


n QQ 
u.yy 


1 77 


iVsi,(4) 


0.16 


2.76 


2.93 


0.16 


3.62 


2.99 


iVsi?(5) 


1.84 


0.16 


0.11 


1.94 


0.56 


0.20 


iVsi3(6) 


1.44 


0.22 


0.16 


1.46 


0.38 


0.18 




0.00 


0.00 


3.31 


0.00 


0.00 


0.00 




15.39 


15.17 


13.46 


15.05 


7.41 


9.89 


A^sd(9) 


5.79 


0.63 


0.44 


6.05 


0.19 


0.82 


\/E'-l[^SD(i)]^ 

y ^ I — \_\. \ /J 


6.44 


5.05 


4.76 


6.31 


3.82 


3.61 


/\ — ^ T Tut /'MO 


6.44 


5.05 


5.80 


6.31 


3.82 


3.61 


yE^=l-6,8,9[^Sfl(i)]2 


17.65 


16.00 


14.28 


17.40 


8.34 


10.56 




17.65 


16.00 


14.66 


17.40 


8.34 


10.56 


1 ri Q n n ("1 1 


/i3 Higgs boson 


iVsi3(l) 


0.00 


0.59 


0.66 


0.01 


0.00 


0.32 




0.00 


0.21 


0.25 


0.00 


0.00 


0.08 




0.00 


0.00 


1.13 


0.00 


0.00 


0.00 


iVsz3(4) 


3.79 


3.43 


3.62 


3.56 


1.55 


4.86 




3.65 


2.51 


2.07 


4.46 


1.54 


1.66 


iVsD(6) 


0.80 


2.13 


1.52 


1.17 


0.38 


1.55 


iVsi3(7) 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 




0.00 


0.00 


9.06 


0.00 


0.00 


0.00 


iVsD(9) 


0.00 


0.77 


0.79 


0.00 


0.00 


0.43 




5.32 


4.80 


4.64 


5.83 


4.76 


5.37 




5.32 


4.80 


4.64 


5.83 


4.76 


5.37 




5.32 


4.86 


10.21 


5.83 


4.76 


5.39 




5.32 


4.86 


10.21 


5.83 


4.76 


5.39 
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Table 3 



Point Number 


1 


2 


3 


4 


5 


6 


Best non-M^VF fusion Nsd 


4.37 {hi) 


3.95 (h2) 


3.62 (/13) 


4.46 (/13) 


4.83 (/ii) 


4.86 (/ig) 


Best WW fusion Nsd 


15.39 (/la) 


15.17 {h2) 


13.46 (/12) 


15.05 (/12) 


16.78 (hi) 


10.08 (/ii) 


Best combined Nsd w.o. 
VFVF-fusion modes 


6.54 (hi) 


5.05 (/i2) 


4.76 (/12) 


6.31 (/12) 


6.69 (/ii) 


5.37 (/13) 


Best combined Ns d with 
M^VF-fusion modes 


17.65 {h2) 


16.00 (/la) 


14.28 (/12) 


17.40 (/12) 


18.07 (/ii) 


10.73 (/ii) 
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